times_preds1 <-read.csv("times_preds1.csv")
par(mar=c(4.1,4.1,2.7,2.1))
plot(times_preds1 $margin~ times_preds1 $at1,type="l",ylim=c(0,1),lwd=5,xlab="Percentage of Member-State Governments Filing Observations",ylab="Probability of Times Coverage of CJEU Decision",col="black",main="",axes=F)
lines(times_preds1 $ci_lb ~ times_preds1 $at1,type="l",lty=2,lwd=2,col="black")
lines(times_preds1 $ci_ub ~ times_preds1 $at1,type="l",lty=2,lwd=2,col="black")
axis(1, c(0,0.1,0.2,0.3,0.4,0.5 ), col = NA, col.ticks = 1)
box(bty="L")
axis(2,c(0.2,0.4,0.6,0.8,1),, col = NA, col.ticks = 1)